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SUMMARY 


A vectorized code, EQUIL, is developed for calculating the equilibrium 
chemistry of a reacting gas mixture on the Control Data STAR-100 computer. The 
code provides species mole fractions, mass fractions, and thermodynamic and 
transport properties of the mixture for given temperature, pressure, and ele- 
mental mass fractions. The code is set up for the e~, H, He, C, 0, N system of 
elements. In all, 24 chemical species are included. 


INTRODUCTION 

This report presents a vectorized code, EQUIL, developed for the Control 
Data STAR-100 computer, which calculates the equilibrium chemistry of a react- 
ing gas mixture. Only gaseous species are considered. The code (see appen- 
dix A) provides species mole fractions, mass fractions, and thermodynamic and 
transport properties of the mixture for given temperature, pressure, and ele- 
mental mass fractions. It can be used as a subprogram to a flow-field code. 

The code is set up for the e“, H, He, C, 0, N system of elements. In all, 

24 species are included. The method given in references 1 and 2 is used in the 
present code for calculating the equilibrium composition. It uses the free- 
energy -minimization technique in which the method of steepest descent is uti- 
lized. Thermodynamic properties, thermal conductivity, and viscosity of each 
input species are evaluated by polynomial expressions. The coefficients for 
these polynomials are prescribed input to the code. The mixture transport 
properties are determined by using the semiempirical relation of Wilke 
(ref. 3) . 

The code is set up for calculating the equilibrium chemistry over 500 mesh 
points. The number of mesh points can be varied by suitably changing the code 
dimensions. The system of elements and species considered in this code can 
also be changed by suitably modifying the input data. 


SYMBOLS 


3] ,a 2 , • • • 
b-| ,b2,b3 
Cl /C 2 

C P 

F 

H 


,a 7 coefficients in approximating 
coefficients in approximating 
coefficients in approximating 
specific heat at constant pressure 
free energy 
enthalpy 


polynomials for 
polynomial for 
polynomial for 


c p . 


U 

k 


H, and F 



k thermal conductivity 

R universal gas constant 

T temperature 

p viscosity 


PROGRAM INPUT 

Most of the data input to the EQUIL code is by punched cards. The 
detailed input information is presented in the following section. 


Chemistry and Transport Model Input 

The elements and species in the mixture are input through NAMELIST THERMO, 
which is described in table I. Thermodynamic properties (refs. 4 to 7) , ther- 
mal conductivity, and viscosity of each input species are represented by 
approximating polynomials. The thermodynamic data are input by cards as 
follows. There are eight cards for each atomic, molecular, or ionic species: 

Card 1 (Format (A6, 4X, 6F5.0)): Field 1 contains the alphanumeric 

identifier for the ith species SYMB(I), which is right justified in 
columns 1 to 6. Fields 2 to 7 (columns 11 to 15, 16 to 20, etc.) 
contain the array AA(I,J) for J = 1 to NE, which specifies the 
number of atoms of each element in the species (see table II) . The 
order must correspond to the order of input of the MHEL array. For 
example, card 1 for species NO + would be 

N0+ -1. 0. 0. 0. 1. 1. 


Cards 2 and 3 (Format (5E14.6)) contain the seven constants a-j to a-j 
used in the polynomials for calculating the specific heat, enthalpy, 
and free energy of each species in the temperature range from 300 K to 
1000 K. These polynomials are 


R 


a-| + a2T + a3T 2 + a^ 3 + 35^ 


H 

RT 


T t 2 a 6 

ai + ao ~ + ao — + a4 — + ac — + — 

1 Z 2 * 3 q 4 D 5 T • 


F 

RT 


a-j (1 - In T) 


a 2 




a 7 


Columns 29 to 80 of card 2 may be used for identification purposes. 
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Cards 4 and 5 (Format (5E14.6)) contain constants a-| to a-j for the 
thermodynamic data in the temperature range from 1000 K to 6000 K. 

Cards 6 and 7 (Format (5E14.6)) contain constants a-j to 87 for the 
thermodynamic data in the temperature range above 6000 K. 

Card 8 (Format (5E14.6)) contains the constants b-| , b 2 , b3, c-| , and 

C2 in the approximating polynomials for viscosity and thermal con- 
ductivity, which are 

U «= b-j + b2T + b3T^ 

k = c-| + c 2 

where T is in K, \i is in lbm/ft-sec, and k is in Btu/ft-sec-°R.* 
Coefficients for the species used in this code are given in appendix B. 

The elements are input in the order e“, H, He, C, O, and N. The order of 
species input is shown in table II. The species composition is described using 
the AA(24,6) array in table II. For example, AA(i,2) defines the number of 
hydrogen atoms in species i, and AA(i,4) defines the number of carbon atoms in 
species i. If species i is a positive ion, then AA(i,l) is -1. 


Temperature, Pressure, and Elemental Mass Fractions Input 

The code is set to calculate the equilibrium chemistry at 500 mesh points. 
This can be varied by changing the dimensions of various variables. The tem- 
perature, pressure, and elemental mass fractions at all the mesh points are 
input to the code through vectors T(500) , P(500) , and CL(500,6). The tempera- 
ture is in kelvins and pressure is in atmospheres. The code requires that at 
least one element of the temperature vector be in each temperature range used 
to calculate the thermodynamic properties in subroutine THERMO. As an example, 
the present code is set to use only two temperature ranges, 1000 K < T < 6000 K 

and T > 6000 K. There is a transition range from 5500 K to 6500 K to assure 

smooth transition from one temperature range to another. The temperature 
vector should have at least one element less than 5500 K, at least one element 

between 5500 K and 6500 K, and at least one element above 6500 K. 

The vector 0.(500,1) is an elemental array and, for the present setup, 
represents the charge balance, which is zero. However, the solution procedure 
does not allow zero values for an elemental array; therefore, the "electron" 
elemental array is set to an arbitrary small number. In this code, 0(500,1) 
is taken as l.E-10 for all cases. No element should have zero mass fraction. 

For cases where the mixture does not include a particular element, the mass 
fraction for that element can be prescribed as an arbitrary small number. 


*1 lbm/ft-sec = 1.488 Pa-sec; 1 Btu/ft-sec-°R = 6226.5 W/m-K. 
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In addition to the aforementioned input quantities, a criterion for the 
convergence of the calculation procedure is required. The convergence crite- 
rion is as follows: 

N 

CRIT = |Xi - Yi| < l.E-6 
i=l 

Here, the values of are the species mole numbers for the current itera- 

tion, and the values of Y^ are the species mole numbers for the previous 
iteration. Thus, when the sum of the absolute values of the changes in the 
mole numbers for all the species is less than 1 .E-6 from one iteration to the 
next, the calculation is terminated. The value of CRIT can be changed if 
necessary. 


CODE STRUCTURE 

Code EQUIL has a main program in which the quantities such as temperature, 
pressure, and elemental mass fractions are prescribed at all the mesh points. 
The value of CRIT is also input here. In addition to the main program, there 
are six subroutines. 

Subroutine READ reads the NAMELIST THERMO and thermodynamic and transport 
data for various species. Subroutine THERMO calculates the thermodynamic prop- 
erties of various species at all the mesh points. Subroutines CHEQ, MINENG , 
and EQSOL use the method of steepest descent to minimize the free energy. An 
initial assumption is made on the mole numbers of various species, and then an 
iterative procedure is followed to find the set of mole numbers of various 
species which minimizes the free energy. 

Knowing the right mole numbers of various species, the subroutine CHEQ 
then calculates the mixture molecular weight, mole fractions, and the enthalpy 
of the mixture. Subroutine TP calculates the mixture specific heat, thermal 
conductivity, viscosity, and Prandtl number. 

Finally, the main program EQUIL converts the mixture viscosity, conduc- 
tivity, and specific heat to SI units. The species mass fractions are also 
obtained here. The printed output includes the mole fractions of various 
species, and the enthalpy, specific heat, viscosity, conductivity, Prandtl 
number, and molecular weight of the mixture. 


LIST OF VARIABLES 

CIS Species mass fraction 

CL Elemental mass fraction 

COND Thermal conductivity of mixture 

CONDI Thermal conductivity of individual species 
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CP Specific heat of mixture 

CPI Specific heat of individual species 

CRIT Convergence criterion 

FORT Free energy of individual species 

HI Enthalpy of individual species 

KTEST A parameter equal to either zero or one 

MOLEF Mole fraction of species 

MW Molecular weight of species 

MWEL Molecular weight of element 

p Pressure 

SH Enthalpy of mixture 

SIG Prandtl number of mixture 

T Temperature 

VIS Viscosity of mixture 

VISI Viscosity of individual species 

WMIX Molecular weight of mixture 

y Mole number of species in current iteration 

YU Mole number of species in previous iteration used to start a new 

iteration for KTEST = 1 

SAMPLE CALCULATIONS 

Two sample calculations are presented here. The first example is for a 
mixture consisting mainly of e", H, He, C, and 0 elements. The input tempera- 
ture, pressure, and elemental mass fractions are listed below: 

T(1 ;100) = 4000. 

T(1 01 jl 00) = 6000. 

T(201 ; 300) = 12570. 

P (1 ; 500) = 6.3549 


5 



CL(1 , 1 ; 500) 
CL (1,3; 500) 
CL (1 ,4;500) 
CL(1 ,5; 500) 


= l.E-10 
= 0.1008 
= 0.4518 
= 0.02443 


CL (1 ,6;500) = 5.E-6 

CL(l f 2;500) = 1. - (Sum of other five elemental mass fractions) 

The mass fraction for the sixth element, nitrogen, was prescribed as a very 
small number. It took 24 iterations for the solution to converge. The actual 
computing time on the Control Data STAR— 100 computer for this example was 
5.3 seconds. The output for this example is given in table III. 

The second example is for air. For this example, the elemental mass 
fractions for H, He, and C were prescribed as small numbers. The input tem- 
perature, pressure, and elemental mass fractions are as follows: 


T(1 ;1 00) = 7000. 

T ( 1 01 ; 1 00) = 6000. 


T(201 ; 300) = 3991.17 


P (1 ; 500) = 1.146 


CL (1,1; 500) = l.E-10 
CL(1 ,2;500) = 5.E-6 
CL (1,3; 500) = 5.E-6 


CL (1 ,4;500) = 5.E-6 


CL(1 ,5; 500) = 0.233 

CL (1,6; 500) = 1. — (Sum of other five elemental mass fractions) 

It took 26 iterations and 6.2 seconds actual computing time for the solution 
to converge. The output for this example is given in table IV. 

In these examples, the input conditions for the first 100 mesh points, 
for the next 100 mesh points, and for the last 300 mesh points were the same; 
consequently, tables III and IV give the results only at mesh points 1, 101, 
and 201 . 
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CONCLUDING REMARKS 


The code presented in appendix A calculates the enthalpy of the mixture 
for given temperature, pressure, and elemental mass fractions. However, in 
most flow-field calculations, it is the temperature which is to be calculated 
for given enthalpy and pressure. To use this code for such calculations, an 
iterative technique can be used in which a temperature is assumed initially. 

For this temperature, the enthalpy of the mixture is obtained from EQUIL. The 
temperature is then changed using the Newton-Raphson technique until the 
enthalpy of the mixture is obtained within desired accuracy to the initially 
known value. This iterative technique can be incorporated in the main program 
EQUIL, and only subroutines THERMO and CHEQ are required to be put in the iter- 
ative loop. 

To save significantly on time for such calculations, a parameter KTEST is 
used. For the initial guess of temperature, KTEST is set equal to zero and for 
subsequent guesses, KTEST is set equal to one. When KTEST =0, the subroutine 
CHEQ starts with a very crude approximation of species mole numbers. For . 
KTEST = 1 , mole numbers calculated in the previous iteration are taken as the 
starting approximation and the solution converges in fewer iterations. 


Langley Research Center 

National Aeronautics and Space Administration 
Hampton, VA 23665 
February 19, 1980 
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PROGRAM LISTING 

The computational program EQUIL is listed in this appendix in Control Data 
STAR-100 FORTRAN language 1.4 (an extension of ANSI FORTRAN for the Control 
Data STAR-100 computer). This program consists of one main program and six 
subroutines. 


PROGRAM EQUIL ( INPUT , OUTPUT , TAPE5= INPUT , TA PF6=0UTPIJT ) 

DIMENSION C I S ( 500 ,24 ) 

COMMON /El/T( 500) ,P( 500) ,WMIX( 500) ,CL( 500, A) ,SH( 500 ) 

C0MM0N/E2/H I ( 500, 24 ) , MOLFF ( 500,24) , CPI ( 500,24) 

C0MM0N/E4/MW ( 24),SYMB(24),AA(24,6),MWFL(6),AAA( 500,24,6) 
C0MM0N/F8/NE,NS 

C0MM0N/F14/VI S ( 500 ) ,COND( 500) ,S IG( 500 ) »CP(500) 
C0MM0N/E15/RR,CRIT,RCRIT,MNS,MMM,MNF,MNFI,MX 
C0MM0N/E16/ENT ( 500,24 ) ,CPI1( 500,24) 

EQUIVALENCE ( C I S ( 1 , 1 ) , ENT ( 1 , 1 ) ) 

REAL MW,MWFL , MOLFF 
RR=1 .987 
CR I T= 1 . E— 6 
RCRIT=.1*CRIT 
M 1 =500 

C T IS IN DEGREE K AMO P IS IN ATM. 

T ( l;100)=4000. 

T( 101 ; 1 00 ) =6000. 

T(201;300)=12570. 

P ( 1 SMI ) =6. 3549 

C THE ELEMENTS ARE IN THE ORDER OF F-,H,HF,C,D AMD N. 

C L ( 1 , 1 :M1 ) = 1.E-10 
CL( 1,2; Ml >=.4229649999 
CL ( 1 , 3 : Ml ) = . 1008 
CL( 1 ,4;M1 ) = .4518 
CL( 1»5:M1 ) = . 02443 
CL ( 1 , 6 : M 1 ) =5 . E— 6 
CALL READ(Ml) 

MNS=M1*NS 
MMM=NE+ 1 
MNE=M1*NE 
MNE 1 =MMM*M 1 
MX=MMM*MMM*M1 
CALL THFRMO ( Ml ) 

KTEST = 0 

CALL CHEQ ( KTEST , Ml ) 

CALL TP(M1 ) 

DO 10 1=1, NS 

C CIS IS THE MASS FRACTION OF SPFCIF I. HI IS IN J/KG. 

CW=MW ( I ) 

CW1=4184./CW 

CIS( 1 ,1 SMI )=M0LEF( 1,1 SM1)*CW/WMIX( 1 SMI ) 

10 HI ( 1,1 :M1 )=HI ( 1, I ;M1 )*CW1 

HERE SH IS IN J/KG, CP IS IN J/KG K, VIS IS IN N SFC/M2 , COND IS 
IN W/M K. 

SH( 1SM1 )=SH( 1;M1 )*4184. 
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CP( l :mi )=CP ( i:mi )*4184./wmi X( l : m i ) 

VIS( l ; Ml )=0.l*viS( l : Ml ) 

COND ( 1 : M 1 )=418.4384*C0ND( 1:M1) 

THE SPECIES ARE IN THE ORDER OF E- , H , H2 ,H + ,HE*HE+ » C , C2 ,C3 , C+ ,C2H , 
C2H2,C3H,C4H,0,02,0+,CO,C02,N,N2,N+,Nn,ND+ 

WR I TE ( 6 , 100 ) 

WRITE (6, 130 ) ( M, MOLEF ( M, 1 ) * MOLEF ( M,2 ) , MOL EF ( M, 3) » MOL EF ( M ,4 ) , MOLFF ( M 
1,5) , M0LEF(M,6) ,M=1,M 1,100) 

WR I TE ( 6 , 1 20 ) 

WRITE! 6,130) ( M , MOL EF ( M , 7 ) , MOLEF ( M,8) ,MOLEF ( M ,9 ) , MOLEF ( M, 10) , 

1M0LEF ( M, 11 ) , MOLEF (M, 12 ) ,M=1,M1 , 100 ) 

WR I TF ( 6 , 1 40 ) 

WRITE (6, 130 ) ( M, MOLEF ( M, 13) , MOLEF! M, 14) , MOL F F ( M , 1 5 ) , MOLF F ( M , 1 6 ) , 
1M0LEF(M,17) , MOLFF (M, 18) ,M=1,M 1,100) 

WR I TE ( 6 , 160 ) 

WRITE! 6,130) ( M , MOL EF ( M , 19 ) , MOLEF ( M ,20 ) , MOLFF ( M , 2 1 ) ,MOLEF(M,22 ) , 
1M0LEF ( M , 23 ) ,M0LEF(M,24) ,M=1,M1, 100) 

WRITE (6, 150) 

WRITE (6, 130) (M,SH(M) ,CP(M) , VIS! M) , COND ( M ) , S I G ( M ) , WM IX ( M ) , M= 1 , Ml , 
1100) 

100 FORMAT ( / ,4X, »M‘ , 1 IX, 'XE' , 13X, ' XH' , 13X, 1 XH2 ' , 1 2X , ' XH+' , 1 2X , • XHE ' , 

1 11X , ' XHE+ 1 ,/ ) 

120 FORMAT! /,4X, 1 M 1 , 1 IX , • XC • , 1 2X , ' XC2 ' , 1 2X , ' XC3 ' , 1 2 X , • XC+ • , 12X, ' XC2H 1 , 
111X, ' XC2H2 ' ,/ ) 

130 F0RMAT(2X,I4,2X,6E15.5> 

140 FORMAT! / ,4X, ' M 1 , 1 1 X , ' XC3H 1 ,11X, ' XC4H ' , 1 1 X , ' XO • , 1 3X , ' X02* ,12X , 'X0+' 
111X, 'XCO' ,/ ) 

160 FORMAT! / ,4X, ' M ' , 11X, 'XC02' , 1 1 X , • XN ' , 1 2X , ' XN2 ' , 1 2 X , 'XN+* ,12X, ' XMO« , 
112X, 'XNO+' ,/ ) 

150 FORMAT! / ,4X, 'M',HX,'SH',12X,'CP',13X,'VIS , ,12X, 'COND' ,11X, ' SIG' , 
111X, 1 WM I X ' ,/) 

STOP 

END 

SUBROUTINE READ! Ml ) 

COMMON/ F3/ A I ( 24 , 3 > , BI ( 24, 3 ) , C I ( 24 , 3 ) , D I ( 24 , 3 ) , F I ( 24 , 3 ) , F I ( 24 , 3 ) , 
1GI ( 24,3) 

C0MM0N/E4/MW! 24 ) , SYMB ( 24 ) , A A ( 24 , 6 ) , MWFL ( 6 ) , AAA ( 500, 24,6 ) 
C0MM0N/E8/NE,NS 

COMMON/ E9/XMA! 24) ,XMR( 24) ,XMC( 24) , XK A ( 24 ) , XKB ( 24 ) 

DIMENSION I LW ( 9 ) , I LWP ( 6 ) 

REAL MW , MWEL 

DATA ILW/5H SI,5H F , 5H 0,5H N,5H C,5H HE, 

15H H,5H E-,5H / 

NAME LI ST/THERMO/ MW , MWEL , ME , MS 

C INPUT PROBLEM NAMELISTS 

20 READ(5, THERMO) 
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55 


60 


WRITE(6, THERMO) 

DO 60 1=1, NS 

REAO( 5 ,901 ) SYMB(I),(AA(I,J),J=1,MF) 

DO 55 J = 1,3 

dca^! I ’« 02 * AI ( 1 » J) » BI < 1 ♦ J) »CI ( 1 t J ) ,01 ( I , J ) ,EI ( I , J ) ,FI ( I , J) ,GI ( I , J ) 

REA0( 5,902 )XMA( I ) ,XMB( I ) ,XMC( I),XKA(I),XKB(I) 

CONTINUE 
C THERMOCHEMICAL PROPERTIES 
WRITE(6,919) 

00 120 1=1, NS 

WRITE (6, 920) SYMB( I ) , ( A I (I , J ) ,8 I ( I , J ) ,C I ( I , J ) , D I ( I , J ) , FI ( I , J ) , 

• FI ( I , J ) , G I ( I,J),J= 1,3) 

120 CONTINUE 

122 WRITE(6,921 ) 

00 130 1=1, NS 

WR 1 TE ( 6 ,922 ) SYMB( I ) ,XMA( I ) ,XMB( I ) ,XMC( I) ,XKA ( I ) , XKB( I ) 

130 CONTINIJF 

C SPECIES/ELEMENTAL COMPOSITION MATRIX 

00 135 1=1, NE 

I LWP ( I ) = I L W ( 1 ) 

1 F ( MWEL ( I). IT. 27.) I LWP ( I ) = I l_W ( ? ) 

IE(MWEL( I >.LT.17. ) I LWP ( I ) = I LW ( 3 ) 

I E ( MWEL (I).LT.15.) I LWP ( I ) = I LW( 4 ) 

I F ( MWEL ( I ) . L T. 1 3 . ) I L W P ( I ) = 1 LW ( 5 ) 

IF(MWEL(I).LT.4.1) ILWP(I)=ILW(6) 

I F ( MWEL ( I ) . LT. 1 . 1 ) I LWP (I ) = I LW ( 7 ) 

IF ( MWEL ( D.LT.0.1) I LWP ( I ) = I LW ( « ) 

I F ( MWEL ( I ) . EO. 0.0 ) I LWP ( I ) = I LW ( 9 ) 

CONTINUE 

WR I TE ( 6 , 923 ) ILWP 
00 140 1=1, NS 
00 137 J= 1 , NE 

IF( AA( I , J ) .EQ.O. ) AA ( I , J ) =0 . 

CONTINUE 

WRITE (6, 924) S YMB ( I ) ,(AA( I ,J) ,J=1,NE) 

CONTINUE 

FORMAT ( A6,4X,6F5.0 ) 

FORMAT { 5E14.6 ) 

F0RMAT(1H1,34X, 'THERMOPHYSICAL PROPERTIES - CURVF FIT COEFFICIENTS 
. ' //46X , ' ( 1 ) THFRMOOYNAMIC PROP ERT I ES ' //4X , 'SPFCIES',11X,'A',14X, 

. ' B ' , 14X , ' C ' , 1 4X , • 0 ' , 14X , • E ' , 1 4X , ' F ' , 1 4X , • G ' / ) 

FORMAT ( 1H0 , 5X , A6 , 3X ,7 El 5. 6, ' T= 300K ' / 1 5X , 7 FI 5. 6 , ' T= 1000K • / 1 5X , 

FORMAT (1H0/4X, 'SPECIES ',30X, 'VISCOSITY' , 2BX , 1H*, 1 2X, 'CONDUCTIVITY* 

• / ) 

FORMAT ( 1H ,4X, A6, AX, 3E 20.6, 3X, 1H«, F 16.6 f F20.6 ) 


135 


137 

140 

901 

902 
919 


920 


921 


922 
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923 FORMAT (1H0// 49 X» 'ELEMENTAL PARTICLES TABLE'/42X,» SPECIES '*6A5) 

924 FORMAT ( 1H ,43X , A6 , IX , 6F5 .0 ) 

NSNEM1=NS*NE*M1 

AAA ( 1 ,1 ,1 JNSNEM1 ),=0. 

DO 110 1=1, Ml 
DO 110 J = 1 ,NS 
DO 110 K=1 , NE 
110 AAA( I ,J»K)=AA(J,K) 


c 


RETURN 

END 

SUBROUTINE THERMO (Ml) 

FREE ENERGY, ENTHALPY, AND SPECIFIC I 
DIMENSION BG ( 500 ) 

COMM0N/El/T(500) ,P(500) ,WMIX( 500) 
COMMON/ E 2/HI ( 500 , 24 ) , MOLEF ( 500,24 
COMMON/ E3/ A I (24,3) ,BI(24,3)»CI(24 
1G I (24,3) 

COMMON/ E5/ FORT (500,24) ,Y( 500,24) , 


FAT BY APPROXIMATING POLYNOMIALS 

CL ( 500 , 6 ) , SH ( 500 ) 

,CPI ( 500,24) 

3) ,DI ( 24,3) ,EI (24,3) ,FI (24,3) , 
(500,24) , YB AR ( 500) 


C0MM0N/E8/NE,NS 

COMMON/ El 2/A9( 500) , A 10 ( 500) ,A13 


( 500) ,AA1 ( 500 ) ,AA2( 500) ,TV1( 500) 


t 


1TV2 ( 500 ) 

C0MMDN/E15/RR,CRIT,RCRIT,MNS,MMM,MNE,MNE1 ,MX 
COMMON/E 16/ENT ( 500,24) , CPU ( 500,24) 

DESCRIPTOR DAI, DA9 
BIT BG 
REAL MOLEF 

CDEFF I Cl NTS ARE INPUT FOR THREE TEMPERATURF RANGES. (1) 300K TO 
1000K , (2) 1000K TO 6000K, AND (3) 6000K TO 15000K. K AND L 

DENOTES THF SET OF COEFFICIENTS THAT ARE BEING USED. COMBINE TO 
ASSURE SMOOTH TRANSITION BETWEEN EACH OF THE THREE TEMPERATURE 
INTERVALS. T IS GENERALLY GREATER THAN 6500K. 

T VECTOR SHOULD CONTAIN ATLEAST ONE FLEMENT LESS THAN 5500K,0MF 
ELEMENT BETWEEN 55005 AND 6500K, AND ONE ELEMENT ABOVE 6500K. 

BG( 1 :M1 )=T( 1 JMl ) .IE. 6500. 

ASSIGN DA9, A9( 1 ;M1 ) 

0A9=Q8VCMPRS ( T ( 1 J M 1 ) , BG ( 1 : Ml ) 5 0A9 ) 

L 1 =08SLEN ( DA9 ) 

BG( 1 JLl )=A9( 1 ;L1 ) .LE.5500. 

ASSIGN DA1,A10( 1 :M1 ) 

DA1=08VCMPRS ( A9( 1 J LI ) ,BG( 1 ;L1 ) 5 DAI ) 

L2=Q8SLEN ( DAI ) 

AA1( 1:L2)=DA1*0A1 
AA2(l:L2)=DAl*AAl(l:L2) 

A13( 1;L2)=AA1 ( 1;L2)#AA1( 1 ; L 2 ) 

TV2( 1 ;L2 )=VALOG( DAI :TV2( 1 :L2) ) 

TV2( 15L2)=1.-TV2( l;L2) 
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no 10 i = 1 , ns 

11*L2)+L* U ,2 )1 BI * 1 t2, * nA l + CI < I»2)*AA1 ( 1 :L2)+ni( I,2)*AA2( 

FORK 1,1 ;L2 )=AI ( I ,2)*TV2( 1;L2)-.5*BI ( I , 2 ) #nA 1 -C I ( I *2 ) *AA1 ( 1 $L2 )/A. 
1 HI ( I ,2 ) *AA2 ( 1 :L2)/12.-EI ( 1 »2)*0.05*A13( 1 1 L 2 ) + FI ( I ,2 ) /0A1-GI (1,2) 
ENT ( 1 , I ; L2 ) =RR*DA 1* ( AI(I,2)+.5*BI(I,2)*DAl+CI(I,2)*AAl(l;L2)/3.+ 
101(1,2 > =:' A A 2 (l:L2)/4.+FI(I,2)*A13(l:L2)/5.+FI(I,2)/nAl) 

CONTINUE 

BG ( 1 ; L 1 ) =A9 ( 1 : l. 1 ) .GT.5500. 

0A1=08VCMPRS(A9 ( 1 ; L l ) , BG ( 1 ; L 1 ) ; OA 1 ) 

L3=08SLEN( DAI ) 

AA1 ( 1 :L3)=DA1*DA1 
AA2( 1 :L3)=DA1*AA1( l;L3) 

A 1 3 ( 1 ; L 3 ) = A A 1 ( 1 :L3)*AA1 ( 1 JL3) 

TV2(l:L3)=VALOG(DAl;TV2( 1 ; L3 ) ) 

TV2(1 :L3)=1.-TV2( 1 : L3 ) 
no 20 1=1, NS 

TV1(1;L3)=RR*( (6.5-.001*0Al)*(AI( I , 2 ) +R I ( I , 2 ) *DA1+C I ( I,2)*AA1(1;L3 

1 ) +0I ( I ,2)*AA2(1!L3)+EI( I,2)*A13< lSL3>>+( .001*DA1-5.5)*(AI( I,3)+BI( 
21 ,3>*0A1+CI(I,3)*AA1(1;L3)+DI( 1 , 3 ) *AA2 ( l ; L3 ) + FI ( 1 , 3 ) *A1 3 ( l ; |_3 ) ) ) 

CPI1(1,I:L1)=Q8VMERG(TV1(1;L3),CPI(1,I;L2),BG(1:L1);CPI1(1,I;L1)) 

TVl(l?L3)=(6.5-.001*DAl)*(AI(I,2)*TV2(l:L3)-.5*BI(r,2)*DAl-Cl(I,2) 

1 -.-A A 1 ( 1 ;L3)/6.-DI ( I , 2 ) *AA2 ( l!L3)/12» — FI ( 1 »2)*.05*A13( 1 ; L3 ) +F I ( I ,2) / 

2DA1-GI ( I , 2 ) ) + ( .001*0Al-5. 5 )*( A I ( I ,3 )*TV2 ( 1 ; L3 )-0. 5*BI ( I ,3 )*0A1-C I ( 

31 , 3 ) * A A 1 ( 1 ;L3 )/6.-DI ( I ,3 >*AA2 ( 1:L3)/12.-FI(I ,3 )*.05*A13( 1 ;L3 ) + 

AFI ( I ,3 )/DAl-GI ( I ,3) ) 

HI < 1 ,1 :H )=08VMERG(TV1( 1 ;L3) , FORT ( 1 , 1 : 12 ) , BG ( 1 ; L 1 ) ; H I ( 1 , I ; L 1 ) ) 

TV1 ( 1 :L3 )=RR*DA1* ((6.5-.001*DAl)*(AI(I,2)+.5*RI( 1,2 )*0A 1+C I (1,2)* 

1 A A 1 ( l:L3)/3. + ni ( I , 2 ) *A A2 ( l;L3)/4. + EI ( I , 2 ) *A1 3 ( l:L3)/5.+FI( I , 2 ) /DAI 

2 ) + ( ,001*DA1-5.5)*(AI ( I , 3 ) + . 5* B I ( I , 3 ) *DA 1+C I ( 1 , 3 )*AA1 ( 1 ;L3 J/3.+0I ( I 
3,3)*AA2(l;L3)/4.+EI(I,3)*A13(l:L3)/5.+FI(I,3)/DAl)) 

^MOLEF ( 1 , I ;L1 )=08VMERG( TV1(1;L3),FNT(1,I:L2),BG(1:L1);M0LEF(1,I;L1) 

CONTINUE 

BG ( 1 : Ml )=T ( 1 ; Ml ) .GT.6500. 

DA1=Q8VCMPRS(T( 1 ? Ml ),BG( 1 ; Ml ) ; DA 1 ) 

L2=Q8SLEN ( DAI ) 

AA1 ( 1 ; L2 )=DA1*DA1 
AA2( 1:L2)=DA1*AA1( l;L2) 

A 1 3 ( 1 ;L2 ) = A A 1 ( 1 :L2)*AA1( 1:L2) 

TV2( 1 ; L2 )=VALOG( DAl;TV2( 1;L2) ) 

TV2( l;L2)=l.-TV2( 1;L2) 
no 30 1 = 1 , NS 

TV 1 ( 1 ; L2 ) =RR* ( A I ( I , 3 ) + B I ( I » 3 ) *DA1+C I ( 1 , 3 ) * A A 1 ( 1 ; L2 ) +0 1(I,3)*AA2(1* 
1L2)+EI ( I , 3 )*A 13 ( 1;L2) ) 

CPI(1,I;M1)=08VMERG(TV1(1;L2),CPI1(1,I:l1),BG(1:M1);CPI(1,I;mi ) ) 
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TV1(1:L2)=AI(I,3)*TV2( 1JL2)-.5*BI ( I * 3 ) *DA 1-C I ( 1 , 3 ) *AA1 ( 1 ; L2 ) /6.-DI 
1 ( I ,3)*AA2( 1 :L2 )/12.-FI ( 1 ,3 ) *0 .05*A 13 ( 1 ; L2 ) +F I ( I ,3)/DAl-GI ( 1,3) 

FORK 1 , I ;M1 )=Q8VMERG(TV1< l:L2> ,HI < 1, I : LI ) ,BG< 1 :M1 ) : FORT ( 1 , 1 JMl ) ) 

TV i ( 1JL2 )=RR*DA1*( AI ( I ,3)+.5*BI ( I,3)*DA1+CI< 1 ,3)*AA1 ( 1 ;L2 J/3.+DI < I 

1,3)#AA2( 1 !L2 J/4.+EI ( I »3)*A13( 1 5L2 J/5.+FI ( I »3)/DAl ) 

ENT ( 1 , 1 5 Ml )=Q8VMERG( TV1 ( 1 ; L2 ) , MOLEF ( 1 , 1 5 L 1 ) »BG ( 1 ;M1 ) ; ENT ( 1 » I ; Ml ) ) 


30 CONTINUE 

TV 1 ( 1 5 Ml ) = VAL OG ( P ( 1 : Ml ) S TV 1 ( 1 ; M 1 ) ) 

DO 50 1=1, NS 

50 FORT ( 1 , I :M1 )=FORT( 1, I :M1 )+TVl ( 1:M1 ) 

HI(1, 1 ;NS*M1 ) =ENT ( 1 ,1 :NS*M1 ) 

RETURN 

END 

SUBROUTINE CH FO ( KT EST ,M1 ) 

C CHEMICAL EQUILIBRIUM OF MULTIPHASE SYSTEMS BASED ON THE PRINCIPLE* 

C OF MINIMIZATION OF THE FREE ENERGY OF THE MIXTURE 

C THE CONDENSED SPECIES OPTION IS NOT CURRENTLY IMPLEMENTED. 
COMMON/El/TI 500) ,P( 500) ,WMI X ( 500 ) ,CL ( 500,6) ,SH( 500) 

COMMON/ E2 /HI (500,24) , MOLEF ( 500,24 ) , CPI ( 500,24 ) 

COMMON/ E 4/ MW (24)»SYMB(24),AA(24,6),MWEL(6),AAA(500,24*6) 

COMMON/ E5/ FORT ( 500,24),Y(500,24),X(500,24),YBAR(500) 

C0MM0N/E6/YI 1 ( 500,24) 

C0MM0N/E8/NE,NS 
COMMON/ E10/SK IP ( 500) 

C0MM0N/E11/C0NV( 500) ,XLAMBD( 500) ,OELT ( 500,24) ,F( 500,24) ,DEBAR( 500) 
1 ,HALL( 500 ) ,DFDL ( 500) 

COMMON/E 1 2/ A9( 500 ) , A10( 500) ,A13(500),AA1( 500 ) , AA2 ( 500 ) , TV1 ( 500 ) , 
1TV2 ( 500) 

C0MMDN/E15/RR,CRIT,RCRIT,MNS,MMM,MNE, MNE1 ,MX 
DIMENSION BIG ( 500) ,WMIX1 ( 500 ) 

EQUIVALENCE ( WMI X 1 ( 1 ) , A9 ( 1 ) ) 

BIT BIG 

REAL MW, MWEL, MOLEF 


NT = 0 

SKIP! 1 :M1 )=0. 

C IF KTEST EQ 1 USE MOLE NUMBERS COMPUTED PREVIOUSLLY FOR THIS 
C STATION AS INITIAL GIJFSS. 

C OTHERWISE , ESTIMATE SPECIES MOLE NUMBERS FROM FLFMENT 
C MASS FRACTIONS. 

I F ( KTEST . EO. 1 )G0 TO 48 

C STARTING ASSUMPTION - ATOMS ONLY, NO COMPOUNDS 
Y ( 1 , 1 ; MNS ) = 1 . E-7 
Y ( 1 , 1 : M 1 I =CL ( 1 , 1 !M1 >/MWEL< 1 ) 

Y( 1,2 5 Ml )=CI.( 1,2:M1)/MWEL(2) 

Y(1,5;M1)=CL(1,3;M1) /MWEL ( 3 ) 

Y( 1,7;M1 )=CL( 1 ,4 ; Ml ) /MWEL ( 4 ) 
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Y(1»15:M1)=CL(1,5;M1) /MWEL ( 5 ) 

Y ( 1 »20 ;M1 ) =CL ( 1 ,6 JMl ) /MWEL ( 6) 

GO TO 50 
48 CONTINUE 

DO 42 I = 1 , NS 

42 Y(1,I;M1)=YI1(1,I;m1) 

50 CONTINUE 

C FREE ENERGY MINIMIZATION BY STEEPEST DFSCFMT 
60 CONTINUE 
NT=NT+ I 

CALL MINENG(NS,NE,M1 ) 

C LAMBDA AND DIRECTIONAL DERIVATIVE (DFDL), AND CONVERGENCE TEST 
XLAMBD ( 1 ; M 1 )=1 . 

DELT( 1 , 1 :MNS )=X ( I , 1 ;MMS)-Y( 1 , I :MNS) 

DO 100 M= 1 t M 1 

IF(SKIP(M).E0.1. )G0 TO 105 
IF(NT.LE.8)G0 TO 107 

I F(NT .GE. 17 .AND. NT. LE. 22 )G0 TO 107 
DO 101 1 = 1 , NS 

IF(DELT(M,I J.GE.O. )G0 TO 101 
I F ( Y 1 1 ( M , I ) • L T . 1 . E— 7 ) GO TO 102 
XL AM=-Y (M, I )/DELT(M, I ) 

IF ( XLAM.GE. XLAMBD! M) )G0 TO 101 
XLAMBD(M)=0. 99999 99 #XLAM 
GO TO 101 

102 DELT ( M , I ) =0 . 

101 CONTINUE 

GO TO 100 
107 DO 108 1 = 1 * NS 

I F ( DELT { M, I ) .GE. 0. ) GO TO 103 
XL AM=-Y ( M , I ) /DELT ( M» I ) 

I F ( XL AM.GE. XLAMBD ( M ) )G0 TO 103 
XLAMBD ( M ) =0.999 999 9* XL AM 

103 CONTINUE 
GO TO 100 

1-05 DO 106 J = 1,NS 
106 DELT ( M , J ) =0. 

100 CONTINUE 

C DERIVATIVE FOR GASEOUS SPECIES. 

F(1 t 1:MNS)=VABS(DELT(1,1;MNS):F(1,1;mnS) ) 

CONV ( 1 : M 1 ) =0 . 

DEBAR ( 1 :M1 ) =0. 

DO 110 1=1 » NS 

CONV( 1 :M1 ) =F ( 1 t I : Ml ) +CONV ( 1 : Ml ) 

110 DEBAR! 1 :M1 )=DELT(1 t I 5M1 )+DFBAR( 1 ;M1 ) 

NTR I ES=0 
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120 CONTINUE 

HALL ( 1 5 Ml ) =1 . / ( YBAR( lJMl ) + XLAMRD( 1 5 Ml ) *DFBAR ( 1 : Ml ) ) 

NTRIES=NTRI ES+1 
00 130 1=1, NS 

130 f(i,i;mi)=(Y(1,i;mi )+xlambo( l: mi )*delt( l , i ; mi ) )*halli i ;mi ) 

F ( 1 , 1 ; MNS ) =VALOG ( F ( 1 , 1 ; MNS ) 5 F ( 1 » 1 : MNS ) ) 

OFDL < l;Ml ) =0. 

00 140 1=1, NS 

F( 1,1 ;M1)=0ELT( 1,1 ;M1 )*( FORK 1,I;M1)+F(1,I:M1)) 

140 DFOL ( 1 ? Ml ) =0FDL ( 1 5 Ml ) +F ( 1 , I 5 M 1 ) 

C IF OFOL < 0, WE ARE GOING THE RIGHT WAY ON FREE ENERGY SURFACE. 

C IF NOT, REDUCE LAMBDA AND TRY AGAIN... 

big(i;mi)=skipu;mi).eq.i. 

OFDL ( 1 ; M 1 )=QBVCTRL( l.E-10,BIG( 1 S Ml ) ; DFDL ( 1 ; Ml > ) 

B I G ( 1 : Ml )=OFDL( 1 ; Ml ).GE. l.E-9 
HALL( 1 5M1 )=.75*XLAMB0( 1:M1 ) 

X L AMBD ( 1 ;M1 ) =08VCTRL ( HALL ( 1 ; Ml ) , R I G ( 1 : Ml ) : XLAMBD ( 1:M1 ) ) 

1 I =08SGE ( DFOL ( 1 5 Ml ) , l.E-9) 

I F ( II.EQ.M1 ) GO TO 200 
IF(NTRIES.GT.16)G0 TO 600 
GO Tn 120 

200 CONTINUE 
C NEW MOLE FRACTIONS 

HALL ( 1 I Ml ) =VABS ( DFDL ( lJMl > ;HALL ( 15 Ml ) ) 

HIG( 1 :M1 )=HALL ( 1 :M1 ) .LT. l.E-9 

CONV ( 1:M1)=Q8VCTRL(RCRIT»BIG( 1 JMl) :COMV( 1 :M1 ) ) 

R I G ( 1 :M1 )=CONV( 1 :M1 ).LE.CRIT 

XLAMBDdJMl )=Q8VCTRL (0. ,BIG( 1 JMl ) J XLAMBD ( 1 :M1 ) ) 

SKI P ( 1 ;M1 ) =QR VCTR L I 1.,BIG(1;M1):SKIP(1:M1)) 

1 1 =Q8SGE (CONV ( 1:M1),CRIT) 

I F ( II.EQ.M1 )G0 TO 600 
DO 220 1=1, NS 

220 y(i,i;md=Y(1,i jmi )+xlambd( i:mi )*delt< l, I ;mi ) 

I F ( NT . LT . 50 ) GO TO 500 
WRITE(6,231 ) 

231 FORMAT! /IX, 'NO. OF ITERATIONS EXCEED 50'/) 

DO 350 M=1 , Ml 

I F ( SKIP(M) .EQ.l. )G0 TO 350 

WRITE (6,300 )M,P( M), TIM) ,DFDL( M) , CONV I M) , XLAMBD! M) 

350 CONTINUE 

300 FORMAT ( IX , 1 4, 5X , • P= • , F 1 2 . 5 , 5X , ' T= ' , F 1 3. 5, 5X , ' DFDL 

1E12.5,3X,'C0NV=',E12.5,3X,'XLAMRD= , ,F11.5) 

GO TO 600 
500 GO TO 60 
600 CONTINUE 

WRITE I 6, 70) NT 
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70 FORMAT! /, 10X , 'NO. OF I TER AT I 0NS = ' , 1 4 , / ) 

VI l ( l , l jmns )=yi l , l : mns ) 

SH ( 1 ; Ml ) =0 . 

WMI XI ( 1 ; Ml ) =0. 


TV2 ( 1 5 M 1 ) =0 . 

00 10 I =1 , NS 
CW=MW ( I ) 

WMI XI ( 1 ;mi ) =WM I x 1 ( 1 JMl )+YI 1( 1, I :M1 > 

10 TV2I 15M1 )=TV2( 1 : Ml )+Y I l ( i,kmi)*cw 
C MOLECULAR WEIGHT OF EQUILIBRIUM MIXTURE 
WMIX I 1 (Ml )=TV2( 1 JMl )/WMIXl( 1 JMl ) 

MIXTURE IN CAL/MOL. DIVIDE BY MIXTURE MOL. 

1 CAL/GM = 4184 J/KG. 

DO 20 I = 1 , NS 


20 


C 

C 


SHI 1 JMl ) =SH ( 1 JMl )+HI ( 1 1 I JMl )#YI 1 ( 1 , I J Ml ) 

MOLFF I 1 , 1 (Ml ) =YI 1 1 1 , I (Ml ) /WMI XI ( 1 J Ml ) 

SHI 1 JMl ) =SH I 1 JMl ) /TV2 1 1 JMl ) 

RETURN 

END 

SUBROUTINE M I NENG I NS ,NE » M 1 ) 

N"s N ^ E Tr E 4 L ^ 4 srt N 4 s T y°sT p E i: T ,N FREE - ENERr ’ v sr « e - 

lxiBAMMO l ? R E H 4 io??;7t’ 6, ’ FI500 ’ 24 ’’ DELT ' 500,24l ’ BI1 " 1500 ’’ 

C0MM0N/E4/MW(24)»SYMB(24) t AA(24»6),MWEL(6)»AAA(500 t 24,6) 

COMMON /E5/ FORT I 500,24) ,Y< 500,24) ,X ( 500,24) ,YBAR( 500) 
COMMON/ E6/YI 1 ( 500,24 ) 

COMMON/ E7/ A 1 ( 500,7,7) ,BB1 1 500,7) 

C0MM0N/E10/SKIPI 500) 

COMMON/ E 1 5/RR ,CR IT ,RCR I T , MMS , MMM ,MME , MMF 1 ,MX 
COMMON/ E 16 /ENT I 500,24) , CP 1 1 1 500,24) 

EQUIVALENCE I F ( 1 , 1 ) , ENT I 1 , 1 ) ) , I OELT I 1 , 1 ) ,C P 1 1 1 1 , 1 ) ) 


DO 100 1=1, NS 

100 YBARI 1 JMl ) =YBAR 1 1 J Ml ) +Y 1 1 , I J Ml ) 
C SET UP AND SOLVE MATRIX 
BB1 (1,1 JMNE1 ) =0. 


DO 110 J=1 ,NE 


DELTA ( 1 , 1 , J j MNS > = AAA( 1 , 1, J JMNS )*Y I 1 , 1 J MNS ) 
DO 110 1=1, NS 

110 BB11 1 , J J M 1 )=DELTA( 1,I,JJM1)+BB1( 1 , J J Ml ) 

C (1) FREE ENERGY - GASEOUS SPECIES 
DO 170 1=1, NS 

170 YI1(1,IjM 1)=Y(1,I;M1) /YBARI 1 J Ml ) 

F(1,1jMNS) = VALOG I Y I 1 ( 1 , 1 j MNS ) J F I 1 , 1 J MNS ) ) 


WT. 
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00 180 1=1, NS 

180 F ( 1 , 1 J Ml ) =Y ( 1 , 1 ; Ml )* ( FORT ( 1 » I : Ml ) + F ( 1 » I J Ml ) ) 

C INITIALIZE MATRICES 
A1( i ,1 ,1 ;mx>=o. 
ri ( l, l , l ; mx )=o. 

A 1 { 1 ,1 ,1 ; MNE ) =BB 1 ( l,l:MNE) 

00 270 J= 1 , NE 
DO 270 K= 1 , J 

OELT( 1,1;MNS)=AAA( 1 , 1 , J : MNS ) *DEITA ( 1, 1,K;MNS) 

00 275 1=1, NS 

275 Rl ( 1 , J,K?M1 ) =OELT (1,1 : Ml ) +R 1 ( 1 , J , K J Ml ) 

270 Rl ( 1 ,K,J;M1 )=R1( 1, J,KJM1 ) 

DO 280 J=2 , MMM 
K = J- 1 

280 A1(1,1,J;MNE)=R1(1,1,K;MNE) 

DO 310 J=2 ,MMM 
K = J- 1 

310 A1(1,MMM,J:M1)=A1( i,k,1:mi) 

DO 320 J=1,NE 

DELT ( 1 ,1 :MNS ) =AAA( 1 , 1 , J ; MNS ) *F ( 1,1; MNS) 

BUM ( 1 ;M1 ) =0 . 

DO 330 1=1, NS 

330 BUM ( 1 JMl )=DFLT( 1 , I ;M1 )+BUM{ 1 ;M1 ) 

320 BBH 1 , J:M1 )=BR1 ( 1, J;M1 )+BUM( 1 ;M1 ) 

BUM ( 1 ; M 1 ) =0 . 

DO 340 1=1, NS 

34 0 BUM (1 JMl ) =BIJM ( 1:M1)+F(1,IJM1) 

BB 1 ( 1, MMM; Ml )=BUM( 1;M1) 

CALL EQSOL ( MMM,NE,M1 ) 

C NEW MOLE FRACTIONS (X) 

XYBAR ( 1 ; Ml ) =BB 1 ( 1 , 1 JMl ) 

BB1 ( 1,1;MNE)=BB1(1,2J MNE ) 

DELT ( 1 , 1 ;MNS )=0. 

DO 390 J=1 ,NE 
DO 400 1=1, NS 

400 DELT( 1,1 JMl )=DELT( 1,1 ;M1 )+AAA( 1 , I , J J Ml )*BB1 ( 1 , J J Ml ) 
390 CONTINUE 

DO 410 1=1, NS 

410 DELT ( 1 , I ; Ml ) =DELT ( 1 , I ; Ml ) +XYBAR ( 1 ;M1 ) 

X( 1 ,1;MNS)=DELT( 1, 1JMNS)*Y( 1 , 1;MMS )-F( 1 , 1 ;MNS ) 

RETURN 

FND 

SUBROUTINE EOSDL < MMM, NE , Ml ) 

COMMON/ E7/ A 1 (500,7,7) ,BB1( 500,7) 

DIMENSION U( 500,7 ) ,S( 500) ,TK2( 500),IJT(500) 

DO 200 K=1 ,NE 
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J = K 

S(i:md=o. 

HO 20 I =K , MMM 

20 S ( 1 JM1 )=S( 1 ; Ml >+A1(1 ,i,j;m1)*a 1(1,1, j : mi ) 

S ( 1 ; Ml ) =VSQRT ( S ( 1 ;mi ) :S( 1?M1 ) ) 

IJ( 1 ,KJM1 ) =A 1 c 1 ,K, J5M1 )+VSIGN( S( 1 : Ml ) , A 1 ( 1 , K , J : Ml ) ; |JT ( 1 ; Ml ) ) 

KP 1=K+1 

KM1=M1*(MMM-K) 

U( 1,KP1;KM1)=A1< l,KPlyJ:KMl) 

ai ( i,k, j:mi )=-vsign(S( i;mi ) »ak i,k, j:mi ) :ijt( i:mi ) ) 

TK2( 1 :M1 )=U( 1 ,K:M1 >*S( 1 5 Ml > 

TK2 ( 1 J M 1 ) =V ABS ( TK 2 ( I J Ml ) ; TK2 ( 1 5 Ml ) ) 

JK=K+1 

DO 70 J J = JK , MMM 
UT ( 1 S Ml ) =0 . 

DO 50 1 1 I =K t MMM 

50 IJT ( 1 ; Ml ) =IJT { 1 ;mi )+U( 1 , 1 1 1 :M1 )*ai ( 1 , 1 1 1 , JJ ;M1 ) 

UT ( 1 ; Ml ) =tJT ( 1 ; M 1 ) /TK2 ( 1 5 Ml ) 

DO AO 1 1 =K * MMM 

40 Al ( l , 1 1 , jj ;M1 ) = a l ( l , 1 1 f JJ5M1 ) — tJ < 1 , 1 1 : Ml )*UT( 1 :mi ) 

70 CONTINUE 

UT ( 1JM1 ) =0 . 

DO 80 1 1 =K » MMM 

80 UT( l ; mi )=ut ( 1 ;mi )+u( i» 1 1 :mi )*bbi ( l » 1 1 :mi ) 

UT( 1;M1)=UT( 1 :M1 )/TK2< 15M1 ) 

DO 60 I I =K » MMM 

60 BB 1 ( 1 y ! I ;M1 ) = BB1 ( 1 * 1 1 :M1 )-u( 1 1 1 1 ;M1 )*IJT( 1 : Ml ) 

200 CONTINUE 

BBH 1»7;M1 )=BB1 ( 1,7: Ml ) /Al ( 1 ,7,7; Ml ) 

BB1(1,6:M1)=(BB1(1,6;M1)-A1( 1 » 6 , 7 J Ml ) *BB 1 ( 1,7:M1))/A1(1,6,6;M1) 
BB1( 1 ,5:M1 )=( BB1 ( 1 ,5:M1 ) —A 1 ( 1 ,5,7? Ml )*BB1 ( 1,7: Ml ) — A 1 ( 1 ,5, 6; Ml )* 
lBBl(l t 6:Ml ) ) / A 1 ( 1,5,5:M1) 

BB1(1,4:M1)=(BB1(1»4;M1)-A1(1,4,7:M1)*BB1(1,7:M1)-A1{1»4,6;M1)* 
1BBH ly6; Ml )-Al ( l,4y 5: Ml )*BB1( 1,5: Ml ) ) / A 1 ( 1,4, 4: Ml ) 

BBH 1 »3:M1 )=( BB1 ( 1 ,3: Ml ) —A 1 ( 1 ,3 ,7; Ml )*BB1 ( 1 ,7:M1 ) —A 1 ( 1 ,3 ,6; Ml )* 
1 BB1 ( 1 1 6 ; Ml ) — Al ( 1 » B » 5 J M 1 ) *BB 1 ( 1 y 5 5 Ml ) — A 1 ( 1 y 3 ♦ 4 S Ml ) HtflBl ( 1 y 4 ; M 1 ) ) / 
2A1 ( 1 ,3,3:M1 ) 

BB1 ( 1,2! Ml ) = (BB1 ( ly 2: Ml ) — Al ( l,2y7:Ml )*BB1( 1 f 7:Ml ) —A 1 ( 1 ,2»6:M1 )* 
1BB1( 1 , 6 ; M 1 ) — A 1 ( 1,2,5;M1 )*BB1 ( ly5;Ml)-Al(l,2,4:Ml)#BBl(l,4;Ml)- 
2A1 ( 1 ,2,3!M1 )*BB1 ( l»3iMl ) )/Al( 1,2,2:M1 ) 
BBl(lyl;Ml)=(BBl(l t l;Ml)-Al(lyly7$Ml)«BBl(ly7:Ml)-Al(lyly6;Ml)* 
lBBl'(ly6;Ml)-Al(lyly5;Ml)#BRl(l,5;Ml)-Al(l,l T 4:Ml)*BBl(l,4;Ml)- 
2A 1 ( 1 , 1 , 3 ; Ml ) #BB 1 ( 1 , 3 * Ml ) — A 1 ( 1 , 1 , 2 J Ml ) *BB 1 ( 1 , 2 J M 1 ) ) / A 1 ( 1 , 1 , 1 ; Ml ) 

RETURN 

END 

SUBROUTINE TP(M1) 


18 



APPENDIX A 


DIMENSION B 1 ( 500 ) , V I S I ( 500,24) , CONDI ( 500,24 ) 

COMMON/El/Tt 500) ,P { 500) ,WMIX( 500 ) ,CL ( 500,6) , SH ( 500 ) 

COMMON/ E2/H I ( 500 , 24 ) , MOLEF ( 500 , 24 ) , CP I ( 500,24) 

C0MM0N/E4/MW( 24) ,SYMB(24) ,AA(24,6) ,MWEL(6) , AAA ( 500 ,24,6) 
C0MM0N/E5/F0RK 500,24) ,Y( 500,24) ,X( 500,24) , YBAR( 500 ) 

COMMON/ EB/NE »NS 

COMMON/ E9/XMA( 24) ,XMB( 24) ,XMC( 24) ,XKA{ 24) , XKB ( 24 ) 

COMMON/ E 1 2 / A9 ( 500) , A 10 ( 500) , A 13 ( 500) ,AA1 ( 500 ) ,AA2( 500) ,TV1 ( 500) , 
IT V2 ( 500 ) 

COMMON/ El 4/ VI S( 500) , CONDI 500 ) ,SIG( 500 ) ,CP ( 500 ) 

COMMON/E 16/ENT ( 500,24) ,CPI 1( 500,24) 

EQUIVALENCE ( VI SI ( 1, 1 ) , ENT ( 1 , 1 ) ) , ( CONDI ( 1 , 1 ) ,CP 1 1 ( 1 , 1 ) ) 

BIT B 1 

REAL MOLEF, MW 

C HERE VISI IS THE VISCOSITY OF SPFCIF I. 

DO 200 1=1, NS 

VISI ( 1 , I ;M1 ) = ( XMA( I )+XMB( I )*T( 1?M1 )+XMC( I )*T< 1 JMl )*T( 1 ?M1 ) )/.0672 
B 1 ( 1:M1)=VISI ( 1,1 ;mi) .LE. o. 

Visit 1,1 ;M1 )=Q8VCTRL( l.E-10,Bl ( l:Ml ) ; VISI ( 1 , 1 :M1 ) ) 

200 CONTINUE 

C CP IS SPECIFIC HEAT OF MIXTURE. 

CP ( 1 : Ml ) =0 . 

DO 210 1=1, NS 

CP ( 1 ; Ml ) =CP ( 1 ;mi )+CPI ( 1, I ;mi )*molef( 1 , I : mi ) 

210 CONTINUE 

C HERE CONDI IS THE THERMAL CONDUCTIVITY OF SPECIE I. 

DO 220 1=1, NS 

CONDI ( 1 » I 5 Ml ) = ( XK A { I ) +XKB ( I ) *T( 1 5 Ml ) ) / .0672 
220 CONTINUE 

C WILKE RELATION FOR MIXTURE VISCOSITY AND THERMAL CONDUCTIVITY, 

v I S ( 1;M1 )=0. 

COND ( 1 ;M1 )=0. 

DO 240 1=1, MS 
A1 3 ( l:Ml ) =0 . 

DO 230 J = 1 , N S 
WJI 1 = MW( J ) /MW ( I ) 

DEN0M1=2.82*S0RT( l.+l./WJI 1 ) 

WJI=S0RT(WJI1 ) 

WJI =SQR T ( W J I ) 

AAl(l;Mi)=viSKl,I:Mi)/viSI(l,J:Mi) 

AAl ( 1 ; Ml ) =VSQRT ( AA1 (1:M1):AA1(1;M1)) 

A A 1 ( 1 : Ml ) =1 , + AAl ( 15M1)*WJI 
AAl ( 1 :M1 )=AA1 ( 1 ;M1 )*AA1( 1 ; M 1 ) 

ai3(1:md=a13(1:mi )+moi.FF( 1 , J:mi )*AA 1( l :mi ) /denomi 
230 CONTINUE 

VI St 1 ;M1 ) =MOL E F ( 1 ,I:M 1 ) * visit 1,I:M1)/A13(1;M1)+VIS(1:M1) 

CONDt 1 ; M 1 ) = M 0 L E F ( 1,1 :M1 )*CONDIt 1,I:M1)/A13(1:M1)+C0ND(1;M1) 

240 CONTINUE 
C PRANDTL NUMBER 

SlGt l :mi )=vi S( l : mi )*CP( l ; mi ) /cond( i:mi)/wmix(1:md 

RETURN 

END 
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PROGRAM INPUT 

The input for the present setup of the program EQUIL is given in this 
appendix. 


f. THERMO 
NS=24,NE=6, 

MWEK 1)=5.486E-4,1.008,4.01 
MW( 1 )=5.486E-4, 1.008,2.016, 
4*01 ,4.01., 

37.041,49.04, 

16. ,32. , 
14. ,28.. 

f.END 

E- 1 

0.2500000E+01 0. 
-0.7453750E+03-0. 1173402E+02 
0 . 2 500000 E+0 1 0. 

-0.745 3749E+03-0. 1173402E+02 
0 . 2508 E+0 1 -0.6332E-05 

-0.7450E+03 -0.1208E+02 

0.0 .0500E-07 

H 1 

0.2500000E+01 0. 
0.2547162E+05-0.4601 176E+0O 
0.2500000E+01 0. 

O. 2547 16 2E+05-0. 4601 176E+00 
2 .47 5 164E + 00 7.366387E-05 

2 . 523626E+04 -3 . 749 1 37E-01 


,12.011,16.0,14.0, 

1 . 008 , 

12.001,24.022,36.033,12.01,25.03,26.038, 

,16. ,28. 011, 44. 01, 

, 14. ,30. ,30. , 


0.294E-05 

H2 


. 0R89E-07 
2 


0.3057445E+01 0.2676520E-02- 
-0. 9889 047 E+0 3-0. 2 29970 5E+01 
0.3100190E+01 0.51 11946E-03 
-0.8773804E+03-0. 1962942E+01 
0.3363E+01 0.4656E-03 

-0. 1018E+04 -0.3716E+01 

-0.079E-05 .0791E-07 

H+ -1 1 

0.2500000E+01 0. 

0. 1840334E+06-0. 1 153R62E+01 
0.2500000E+01 0. 

0. 1840334E+06-0. 1 153862E+01 
0 . 2500000 E+01 0. 

0. 1840334E+06-0. 1 153862E+01 
0.0 .0500E-07 

HE 1 

0.2500000E+01 0. 

-0.745 3749 E+03 0 . 9 153488 E+00 
0.2500000E+01 0. 


0 . 0 . 0 . 

GORDON AMD MCBRIDF NASA SP-273 

0 . 0 . 0 . 

GORDON AMD MCBRIDE NASA SP-273 

0.1 364E-08 -0.1094E-12 0.2934F-17 

ESCH ETAL NASA CR-1119R9 
— . 1000R-12 26.000E-05 0.0 

0 . 0 . 0 . 

GORDON AND MCBRIDE NASA SP-273 

0. 0. 0. 

GORDON AMD MCBRIDE NASA SP-273 

— 2.537593 E-08 2. 3R6674F-12 -4 . 55 143 1 E-17 

NICOLET NASA CR-132470 

—.081 1 E— 1 2 2.496E-05 5.1290E-08 

0.5809916F-05 0 . 552 1039 F-08-0 . 1 R 12273E-1 1 
GORDON AND MCBRIDE NASA SP-273 
0.5264421 E-07-0. 3490997 E-10 0.3694534E-14 
GORDON AMD MCBRIDE NASA SP-273 

— 0. 5 1 27 F-07 0.2R02E-11 -0.4905E-16 

ESCH ETAL. NASA CR-1119R9 
-.0BB6E-12 3.211E-05 5.3440E-0R 

0. 0. 0. 

GORDON AND MCBRIDE NASA SP-273 
0. 0. 0. 

GORDON AND MCBRIDE NASA SP-273 
0 • 0 . 0 • 

GORDON AND MCBRIDF NASA SP-273 

-. 1000E— 12 26.000E-05 0.0 

0. 0. 0. 

GORDON AMD MCBRIDE NASA SP-273 

0 . 0 . 0 • 
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—0*74537495+03 0.91534885+00 GORDON AND MCBRIDE NASA SP-273 

0. 25000005+01 0. 0. 0. 0. 

-0 *74 5 3749 E+ 03 0 .91 53488 E+00 GORDON AMD MCBRIDF NASA SP-273 

-1.34515-06 2.3119 19E-08 -4. 7359 flflF-1 3 2.038RF-05 3.24935-08 

HE+ -1 1 

0 • 2 500000E + 0 1 0. 0. 0. 0. 

0.28534265+06 0.16084045+01 GORDON AND MCBRIDE NASA SP-273 

0.25000005+01 0. 0. 0. 0. 

0.2853426E+06 0.16084045+01 GORDON AMD MCBRIDE NASA SP-273 

0.25000005+01 0. 0. 0. 0. 

0.28534265+06 0.16084045+01 GORDON AMD MCBRIDE NASA SP-273 

0.0 .05005-07 -.10005-12 26.0005-05 0.0 

C 1 

0.25328705+01-0.1588764E-03 0.30682085-06-0.26770065-09 0. 87488825-13 
0.8524042E+05 0.46062375+01 GORDON AND MCBRinE NASA SP-273 
0.2581066E+01-0. 1469620E-03 0. 74388085-07-0 . 7948 107E-1 1 0. 58900975-16 
0.85216295+05 0.4312887E+01 GORDON AND MCBRIDF NASA SP-273 

0.2141E+01 0.3219E-03 -0.54985-07 0.36045-11 -0.55645-16 

0.85425+05 0.6874E+01 ESCH 5TAL NASA CR-111989 

1 • 997E-05 .17725-07 -.33785-12 2.506E-05 .74795-08 

C 2 2 

0. 7451 8 14E+0 1-0.10144685-01 0.85879735-05 0.87321105-09-0. 24429795-11 
0.9891198E+05-0. 15846675+02 GORDON AND MCBRIDF NASA SP-273 
0.40435355+01 0 . 2057 36 5E-03 0. 1 0907 57F-06-0 . 36427875- 1 0 0.34127865-14 
0.9970948E+05 0 . 1 2775 1 5E+01 GORDON AMD MCBRIDE NASA SP-273 

0.4026E+01 0.4857E-03 -0.7026E-07 0.4666E-11 -0.1142E-15 

0 . 9787E + 05 0.10905+01 ESCH ETAL NASA CR-111989 

1.931E-05 .1393E-07 -.25755-12 .859E-05 .62335-08 

C3 3 

0.57408465+01-0. 84281 2 3E-02 0. 18620195-04-0.14510525-07 0. 39676975-11 
0.97157525+05-0.23837375+01 GORDON AMD MCBRIDE NASA SP-273 
0.36815365+01 0. 241 6 52 35-02-0. 84348 1 1 F-06 0.14508195-09-0.95697305-14 
0.97413955+05 0.68377805+01 GORDON AND MCBRIDF NASA SP-273 

0.22 1 3E+02 -0.17595-01 0.55655-05 -0.67585-09 0.28255-13 

0.9423E+05 -0.10215+03 ESCH ETAL. NASA CR-111989 

2.019E-05 .1179E-07 -.16555-12 .6305-05 .58045-08 

C+ -1 1 

0.2595384E+0 1-0. 40686645-03 0. 6892366F-06-0. 52664875-09 0.15083375-12 
0.21666285+06 0 . 3895729 E+0 1 GORDON AND MCBRIDF NASA SP-273 
0.2511 827E+0 1-0. 1735978 E-04 0.95042675-08-0. 2218851 5-1 1 0.18621895-15 
0.21667725+06 0 .42 861 29E+0 1 GORDON AMD MCBRIDF NASA SP-273 

0.25285+01 0.4869E-05 -0.70265-08 0.1134E-11 -0.34765-16 

0.2 168E+06 0.41395+01 ESCH ETAL. NASA CR-111989 

0.0 • 0500E-07 -.10005-12 26.0005-05 0.0 

C2H 1 2 

0.2649940E+01 0 . 8 49 19 5 1 E-02-0. 98 1 6537 5-0 5 0 . 6537362 E-08-0. 1 7356275-1 1 


HE 

HF 

HE 

HF 

H5 + 

H 5 + 

H5 + 

H 5 + 

H5 + 

H5 + 

HF + 

C 

c 

c 

c 

c 

c 


C2 
C 2 
C2 
C2 
C2 
C 2 


C3 
C3 
C3 
C3 
C3 
C 3 


C + 

c+ 

c+ 

c+ 

c+ 

c+ 


C2H 


5 

6 
7 
R 

2 

3 

4 

5 

6 

7 

8 
1 
2 

3 

4 

5 

6 

7 

8 
1 
2 

3 

4 

5 

6 

7 

8 
1 
2 

3 

4 

5 

6 

7 

8 

2 

3 

4 

5 

6 

7 

8 

2 
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0.5627575E+05 0.7689860E+01 GORDON ANO MCBRIDE NASA SP-273 
0.4420765E+01 0 .22 1 1930E-02-0.5929494E-06 0.9419577E-10-0.6R52759E-14 
0.55835 44E+05-0. 11588 09 E+01 GOROON AMD MCBRIDE NASA SP-273 

0.5307E+01 0.8966E-03 -0.1378E-06 0.9251E-1I -0.2278F-15 

0. 5809E+05 -0. 5288E+01 ESCH ETAL. NASA CR-11I989 

2.404E-05 .1363E-07 -.2184F-12 1.126E-05 .7439F-08 

C2H2 22 

0.14 10276 E+01 0.1 90 5727 E-0 1-0. 24501 39 F-04 0. 1639087E-07-0.4134544E-1 1 
0.2618820E+05 0. 1 139382E+02 GORDON AND MCBRIDE NASA SP-273 
0.4575108E+01 0.5123835E-02-0.1745235E-05 0 . 2867306E-09-0 . 1795142E-13 
0.2560742E+05-0. 3573794E+01 GOROON AND MCBRIOE NASA SP-273 

0.6789E+01 0. 1503E-02 -0.2295E-06 0.1534E-10 -0.3763E-15 

0.2590E+05 -0. 1 539E+02 ESCH ETAL. NASA CR-111989 

1.396E-05 .0842E-07 -.6939F-12 1.126E-05 .7439E-08 

C3H 1 3 

3. 3446607 E+00 1 .06 87605 E-02-1 .3312 138F-05 1 . 33R9601E-0R-5.6987727E-12 
6.2581906E+04 6 . 00041 84E+00 WAKELYN AND MCLAIN 72657 

3. 8776821 E+00 6 .7242969E-03-2 .6055734F-06 4.4163330E-10-2 .7082704E-14 
6.2564338E+04 3 . 8 265297E+00 WAKELYN AND MCLAIN 72657 

3. 8776821 E+00 6.7242969E-03-2.6055734F-06 4. 4 1 633 30 F- 10-2 .7082704E-14 
6.2564338E+04 3 . 8 265297 E+00 WAKELYN AND MCLAIN 72657 

2.019E-05 . 1 179E-07 -.1655F-12 .630E-05 .5804E-08 

C4H 1 A 

4.96866 10 E+00 1 .7278593 E-02-2 .9943 17 1 F-05 3.2461613E-08-1 .3663978E-11 
7 . 5454605E+04-8 .7699380E-01 WAKELYN AND MCLAIN 72657 

6 .53 1 2534E+00 6 . 506462 1 E-03-2 . 2 51 741 1 E-06 3. 329 57 82 F- 10-1 .7214711F-14 
7. 535041 2E+04-7 . 4467228E+00 WAKELYN AND MCLAIN 72657 214 

6 . 53 1 2534E+00 6 . 506462 1 E-03-2 .2 51741 1 F-06 3.3295782E-10-1.7214711E-14 
7 . 53504 1 2 E+0 4-7 .4467228E+00 WAKELYN AND MCLAIN 72657 214 

2.019E-05 . 1.179E-07 -.1655F-12 .630E-05 .5804F-08 

0 1 

0. 294642 8 E+0 1-0. 16 38 166 E-02 0 .242 1 03 1 F-05-0 . 1602843E-08 0 .3890696E-12 
0.29 14764E+05 0. 2963994E+01 GORDON AND MCRRIDE NASA SP-273 
0.2542059E+01— 0.2755061 E-04—0 .310280 3F- OR 0 . 45 51067 F- 1 1 -0 .4368051 F-l 5 
0. 2923080E+05 0. 4920308E+01 GORDON AND MCBRIDE NASA SP-273 

0 . 2546E+0 1 -0.5952E-04 0.2701E-07 -0.2798F-11 0.93R0F-16 

0.29 1 5F+05 0. 5049E+01 ESCH ETAL. NASA CR-111989 

1 • 5 19E-05 . 1875E-07 -.222RF-12 1 .250E-05 .7092F-0R 

02 2 

0. 362 5598 E + 01 -0. 18782 18 E-02 0.7055454F-05-0 .67635 13 E-08 0 .2 155599E-1 1 
-0. 1047522E+04 0 . 4305277E+0 1 GORDON AND MCRRIDF NASA SP-273 
0.362 1953E + 01 0 .736 1 8 26 E-0 3-0 . 196 5 22 2 E-06 0.36201 55F-10-0. 2894562 F- 14 
-0. 1201982E+04 0 . 36 1 5096E+0 1 GORDON AND MCBRIDE NASA SP-273 

0.372 1 E+0 1 0.4254E-03 -0.2835E-07 0.6050F-12 -0.51R6F-17 

— 0. 1 044E+04 0. 3254E+01 ESCH ETAL. NASA CR-111989 

1.693E-05 . 1 496 E-07 -.2276E-12 1.019E-05 .4901E-0R 


C2H 3 

C2H 4 

C2H 6 

C2H 7 

C2H 8 

C2H2 2 

C2H2 3 

C2H2 4 

C2H2 5 

C2H2 6 

C2H2 7 

C2H28 

C3H 2 

C3H 3 

C3H 4 

C3H 5 

C3H 6 

C3H 7 

C3H 8 

C4H 2 

C4H 3 

C4H 4 

C4H 5 

C4H 6 

C4H 7 

C4H 8 

1 

0 2 

0 3 

0 4 

0 5 

0 6 

0 7 

0 8 

1 

02 2 

02 3 

02 4 

02 5 

02 6 

02 7 

02 8 
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0 + -1 1 

0.24984796+01 0. 1 141 097E-04-0. 29761 396-07 0.32246586-10-0. 12375516-13 0+ 

0.1 879490 E+06 0 .4386435 E+0 1 GORDON ANO MCBRIDE NASA SP-273 0+ 

0. 2 5 06 048 E +01-0. 1446424E-04 0.12446046-07-0.46858476-11 0.65548876-15 0+ 

0. 1879470E+06 0.4347974E+01 GORDON AND MCBRIDE NASA SP-273 0+ 

0. 2944E+0 1 -0.4108E-03 0.9156E-07 -0.5848E-11 0.1190E-15 0+ 

0.1879E+06 0.17506+01 ESCH ETAL. NASA CR-111989 0+ 

0.0 .05006-07 -. 1000E-12 26.0006-05 0.0 0+ 

CO 1 1 

0.3710092E+01-0. 1619096E-02 0.36923596-05-0.20319676-08 0.23953346-12 CO 
-0.1435631E+05 0.2955535E+01 GORDON AND MCBRIDE NASA SP-273 CO 

0.2984069E+01 0 . 1489 1 39E-02-0. 57899686-06 0. 10364576-09-0. 6935355E-14 CO 
— 0. 1424522E + 05 0 .63479 1 5 E+.0 1 GORDON AND MCBRIDE NASA SP-273 CO 

0.3366E+01 0.8027E-03 -0.19686-06 0.19406-10 -0.55496-15 CO 

-0.1434E+05 0.42636+01 ESCH ETAI.. NASA CR-111989 CO 

2.404E-05 . 1 363E-07 -.21846-12 .8596-05 .62336-08 CO 

C02 1 2 


0.24007796+01 0 . 873509 5E-02-0. 66070876-0 5 0.20021866-08 0.63274036-15 C02 

—0.48 377 5 2 E+0 5 0.96951456+01 GORDON AND MCBRIDE NASA SP-273 C02 

0. 44608 04E + 01 0. 3098 17 IE-02-0. 12392576-05 0.22741326-09-0. 15525956-13 C02 

—0 .48961 44E+05-0 .9863598 E+00 GORDON AND MCBRIDE NASA SP-273 C02 

0. 4413E+0 1 0. 3 192E-02 -0.12986-05 0.24156-09 -0.16746-13 C02 

-0.48946+05 -0 . 72 88 E+00 ESCH ETAL. NASA CR-111989 C02 

2 . 404E-05 .13636-07 -.21846-12 .8596-05 .62336-08 'CO 

N 1 

0. 25030716+01-0. 2180018E-04 0. 54205286-07-0 . 5647 560E- 1 0 0.20999046-13 N 
0. 5609890 E+05 0 .41 67576 E+0 1 GORDON AND MCBRIDE NASA SP-273 N 


0.2450268E+01 0 . 1 066 14 5E-03-0. 74653376-07 0.18796526-10-0.10259836-14 N 
0. 5611604E+05 0 . 4448758 E+0 1 GORDON AND MCBRIDE NASA SP-273 N 

0. 2748 E + 0 1 -0. 3909E-03 0.133BE-06 -0.1191E-10 0.33696-15 N 

0 . 5609E+05 0.2872E+01 ESCH ETAL. NASA CR-111989 N 

0.253E-05 . 2206E-07 -.37376-12 1.2816-05 .85936-08 N 

N2 2 

0. 3674826E+0 1-0. 1208 150E-02 0. 23240 10E-05-0. 632 1 755F-09-0. 22 577256-1 2 N2 
-0.1061158E+04 0.23580426+01 GORDON AND MCBRIDE NASA SP-273 N2 

0.28963196+01 0. 1515486E-02-0. 57235276-06 0.99807396-10-0.65223556-14 M2 
-0.90586186+03 0.61615146+01 GORDON AND MCBRIDE NASA SP-273 N2 

0 . 3727E+0 1 0.46846-03 -0.11406-06 0.11546-10 -0.3293F-15 N2 

-0 . 1043E+04 0.12946+01 ESCH ETAL. NASA CR-111989 N2 

0.970E-05 .16 1 3E-07 -.19166-12 .6546-05 .6457F-0R N2 

N+ — 1 1 

0.2727E+01 -0.2820E-03 0.11056-06 -0.15516-10 0.78476-15 N+ 

0 . 2254E + 06 0.3645E+01 ESCH ETAL. NASA CR-11989 N + 

0.2727E+01 -0. 2820E-03 0.11056-06 -0.15516-10 0.78476-15 N+ 

0.2254E+06 0.36456+01 ESCH ETAL. NASA CR-111989 N+ 

0. 2499E+0 1 -0.37256-05 0.11476-07 -0.11026-11 0.30786-16 N+ 

0. 2254E+06 0.49506+01 ESCH ETAL. NASA CR-111989 N+ 

0.0 .05006-07 -.10006-12 26.0006-05 0.0 N+ 

MD 1 1 

0.4045952E+01-0. 3418178 E-02 0.79819196-05-0.61139316-08 0.15919076-11 ND 
0.9745393E+04 0. 2997499E+01 GORDON AND MCBRIDE NASA SP-273 ND 

0.3 1 89000E+0 1 0 . 1338228E-02-0. 5289932 E-06 0.9591933E-10-0.6484793E-14 ND 
0.9828329E+04 0.67458136+01 GORDON AND MCBRIDE NASA SP-273 ND 

0.3845E+01 0.2521 E-03 -0.26586-07 0.21626-11 -0.63816-16 NO 

0.9764E+04 0.32 1 2E+01 ESCH ETAL. NASA CR-111989 NO 

0.970E-05 .16136-07 -.19166-12 .654E-05 .64576-08 NO 

N0+ -1 1 1 

0. 36685066+0 1-0. 11 54458 E-02 0.2 17556 IE-05-0. 48227476-09-0. 2784791 6-1 2 N0+ 

0. 1 180337E+06 0. 31779326+01 GORDON AND MC8RID6 NASA SP-273 N0+ 

0.2888549E+01 0.1521712E-02-0.5753124F-06 0.10051086-09-0.66044296-14 N0+ 

0. 1181924E+06 0.70027206+01 GORDON AND MCBRIDE NASA SP-273 N0+ 

0.2888549E+01 0.1521712F-02-0. 57531246-06 0.10051086-09-0.66044296-14 N0+ 

0.118 1924E+06 0.70027-206+01 GORDON AND MCBRIDE NASA SP-273 N0+ 

0 .9706-05 .16136-07 -.19165-12 .6546-05 .64576-08 N0 + 


2 

3 

4 

5 

6 

7 

8 
1 
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TABLE I.- NAMELIST THERMO 


Variable 

Definition 

NS 

Number of species present in the mixture; NS = 24 

NE 

Number of elements present in the mixture, including electrons; 


NE = 6 

MWEL 

An array of molecular weights of the elements present 

MW 

An array of molecular weights of the species present 













































TABLE III.- RESULTS FOR FIRST EXAMPLE 


> Species (i) 

Mole fractions for mesh point - 

1 

101 

201 

e” 

• 1 0356E-1 0 

.65454E-04 

.5521 6E-01 

H 

• 41 791 E+00 

•85288E+00 

.78588E+00 

H 2 

•441 88E+00 

.1 6587E-01 

.10254E-03 

H + 

. 1 0030E-06 

• 59900E-05 

•33278E-01 

He 

• 81 824E-01 

.531 00E-01 

•49085E-01 

He + 

.n 2 05E- 2 0 

• 89638E-1 5 

• 33561E-06 

C 

• 23667E-02 

• 71 937E-01 

.51 611 E-01 

C 2 

•71 464E-03 

• 1 7929E-02 

• 1 4203E-05 

C 3 

•94738E-03 

• 22540E-04 

.98840E-1 1 

c + 

•65849E-06 

. 60351 E-04 

• 21 840E-01 

c 2 H 

.251 02E-01 

•31818E-03 

.1 0560E-08 

c 2 h 2 

•14824E-01 

• 30029E-05 

• 24556E-1 3 

c 3 h 

.6831 3E-02 

• 34229E-05 

.1 311 6E-1 4 

c 4 h 

• 2621 9E-02 

.11 006E-06 

.66449E-1 8 

0 

•47058E-07 

•67763E-04 

• 28840E-02 

°2 

•67916E-15 

. 11 41 8E-1 0 

• 51 267E-09 

0 + 

•97430E-14 

.1 0728E-08 

.98447E-04 

CO 

• 49752E-02 

• 31 593E-02 

.661 40E-06 

co 2 

•2041 6E-09 

•75944E-09 

.5391 5E-1 4 

N 

• 1 0220E-05 

• 1 8098E-05 

• 1 5684E-05 

n 2 

• 21 336E-05 

.36198E-09 

.57447E-1 4 

N+ 

•72752E-1 3 

.9351 0E-1 1 

• 1 0534E-06 

NO 

.261 04E-10 

.11 051E-09 

•26882E-1 1 

NO + 

•89366E-1 2 

• 1 9795E-1 1 

.90302E-12 




TABLE III.- Concluded 


Thermodynamic and 
transport quantities 
for the mixture 

Mesh point 
1 

Mesh point 
101 

Mesh point 
201 

Enthalpy, J/kg 

.68375E+08 

. 1 7007E+09 

. 28681 E+09 

Specific heat, J/kg-K 

•97937E+04 

. 1 0073E+05 

.12513E+05 

Viscosity, N-sec/m^ 

.67567E-04 

.1111 3E-03 

.1 7887E-03 

Conductivity, W/m-K 

.13222E+01 

. 1 8599E+01 

. 34679E+01 

Prandtl number 

.50052E+00 

.601 92E+00 

.64548E+00 

Mixture molecular weight 

. 32553E+01 

.2111 8E+01 

.1 9523E+01 











TABLE IV.- RESULTS FOR SECOND EXAMPLE 


Species (i) 

Mole fractions for mesh point - 

1 

101 

201 

e" 

• 51 200E-03 

.17932E-03 

.76507E-06 

H 

.10608E-03 

• 1 2973E-03 

• 1 4546E-03 

H 2 

. 63731 E-10 

.31662E-10 

.78930E-08 

H+ 

• 521 21 E-08 

• 1 6868E-08 

• 97069E-1 0 

He 

. 24459E-04 

. 2 9650E-04 

•33909E-04 

He + 

.301 39E-14 

. 1 0362E- 2 0 

• 1 2286E-23 

C 

• 22278E-04 

.54805E-06 

•2221 7E-1 0 

C 2 

.251 86E-11 

. 1 71 95E-1 2 

• 33864E-1 9 

C 3 

.64357E-23 

.51 556E-21 

.13274E-25 

c + 

• 36922E-06 

• 1 9954E-08 

• 25887E-1 3 

c 2 H 

.311 34E-18 

.1 6576E-22 

.74580E-23 

C 2 h 2 

.1 8587E-27 

.11 607E-24 

• 57230E-25 

c 3 h 

. 33961E-23 

. 38821 E-25 

.501 37E-29 

c 4 h 

• 21 375E-30 

. 37541 E-32 

• 20308E-40 

0 

•261 23E+00 

• 31 1 68E+00 

. 25236E+00 

°2 

• 33795E-04 

• 24698E-03 

• 34585E-01 

0 + 

.53787E-04 

• 36374E-05 

.4971 6E-09 

CO 

. 25821 E-04 

.61 045E-04 

.61 969E-04 

co 2 

.81 239E-09 

• 1 0493E-07 

•89826E-06 

N 

• 47790E+00 

• 1 6095E+00 

.1 3023E-02 

n 2 

• 25688E+00 

.51 837E+00 

.66646E+00 

N+ 

• 1 2045E-03 

• 1 6831 E-05 

.11 455E-09 

NO 

.27366E-02 

.81 706E-02 

.45044E-01 

NO+ 

• 34498E-03 

• 1 8364E-03 

• 1 2500E-04 
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TABLE IV.- Concluded 


Thermodynamic and 
transport quantities 
for the mixture 

Mesh point 
1 

Mesh point 
101 

Mesh point 
201 

Enthalpy, J/kg 

.2561 2E+08 

.1 4583E+08 

•73277E+07 

Specific heat, J/kg-K 

.15926E+04 

. 1 3997E+04 

.1 3339E+04 

Viscosity, N-sec/ra 2 

.1921 6E-03 

•15987E-03 

.11 002E-03 

Conductivity, W/m-K 

.39288E+00 

.31 231 E+00 

.21 043E+00 

Prandtl number 

. 77901 E+00 

.71 658E+00 

.69747E+00 

Mixture molecular weight 

.1 81 60E+02 

. 2201 5E+02 

.251 77E+02 
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